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Abstract. We consider a branching random walk initiated by a single 
particle at location in which particles alternately reproduce according 
to the law of a Galton-Watson process and disperse according to the 
law of a driftless random walk on the integers. When the offspring 
distribution has mean 1 the branching process is critical, and therefore 
dies out with probability 1. We prove that if the particle jump distribution 
has mean zero, positive finite variance tf, and finite 4 + £ moment, and if 
the offspring distribution has positive variance and finite third moment 
then the distribution of the rightmost position M reached by a particle of 
the branching random walk satisfies PjM > x] ~ 6rf/{a^x^) asx —> oo. We 
also prove a conditional limit theorem for the distribution of the rightmost 
particle location at time n given that the process survives for n generations. 



1. Introduction 

It has been known since the work of McKean [171 that the time evolution 
of a one-dimensional branching Brownian motion is intimately connected 
with the behavior of solutions of the Fisher-KPP equation 

McKean observed that the cumulative distribution function of the position 
Kf of the rightmost particle at time Hn a one-dimensional branching Brow- 
nian motion obeys equation ([T]). In the simplest case, where particles move 
independently along brownian trajectories and undergo simple binary fis- 
sion following independent, exponentially distributed gestation times, the 
fimction / is given by /(m) = — u, the case originally studied by Fisher 
||12]; more general branching mechanisms lead to more general functional 
forms. In general, when the underlying branching process is supercritical, 
the solution of (H) with Heaviside initial data approaches a traveling wave 
with a positive asymptotic velocity, and thus, in particular, the distribution 
of Rf, when centered at its median mt, converges weakly to the distribution 



Date: December 13, 2012. 

1991 Mathematics Subject Classification. Primary 60J80, secondary 60J15. 
Key words and phrases, branching random walk, critical branching process, nonlinear 
convolution equation, Feynman-Kac formula. 

First author supported by NSF grant DMS - 0805755. 

1 



2 



STEVEN P. LALLEY AND YUAN SHAO 



described by the wave. There is now a considerable literature surrounding 
problems connected with this phenomenon: see, e.g., [8] for the precise as- 
ymptotic behavior of the function 1 1-> mt} [15] for a proof that the traveling 
wave is a mixture of extreme value distributions; and [IJ and [2] for proofs 
that the distribution of the entire point process of particle locations, when 
centered at mt, converges in law as f — > oo. Generalizations of some of these 
results to supercritical branching random walks are given in Q and [Zl- 

When the branching mechanism is critical (that is, when the mean number 
of offspring particles at a particle death is 1), the nature of the process Rt is 
entirely different, because in this case the process must ultimately die out, 
with probability one. Hence, it is more natural in the critical case to ask 
about the distribution of 

(2) M:=maxRf, 

t<oo 

the rightmost point ever reached by a particle of the branching process. 
Interest in the distribution of M stems in part from its relevance to evolu- 
tionary biology, where critical branching Brownian motion has been used 
as a model for the spatial diffusion of alleles with no selective advantage 
or disadvantage: see, for instance, [9], [13J, and the references therein. For 
critical (or subcritical) branching Brownian motion the distribution function 
w{x) - P{M < x} of M satisfies the ordinary differential equation 

(3) ^zv"{x) = -xpiwix)), 

where ipiw) is the probability generating function of the offspring distribu- 
tion of the branching process. This differential equation is similar to that 
satisfied by the traveling wave(s) for the Fisher-KPP equation, but differs 
in the nonlinear term ip{w); this leads to solutions of a very different charac- 
ter, reflecting the qualitative differences between critical and supercritical 
branching. Using the equation (|3|, Fleischman and Sawyer [13J proved that 
if the offspring distribution has mean 1, positive variance, and finite third 
moment then 

(4) P{M >x} ^ as x^ oo. 

a-^x^ 

(In the case of "double-or-nothing" branching, where particles produce 
either or 2 offspring with probability 1 /2, the solution to the differential 
equation with the appropriate boundary conditions has the closed form 
1 - mix) = 6/(1 + xf.) 

Our primary objective in this paper is to show that under suitable hy- 
potheses the asymptotic formula (|4| holds generally for critical, driftless 
branching random walks. For ease of exposition we will limit our study to 
discrete-time branching random walks on the integer lattice Z in which the 
reproduction and dispersal mechanisms are independent, with dispersal 
preceding reproduction. Thus, in each generation, particles 
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(A) jump (independently of one another) to new sites, with jumps fol- 
lowing a mean-zero, finite variance, jump distribution f rw {^xlxez; 

(B) then reproduce as in a simple Galton-Watson process, according to 
a fixed offspring distribution few := {pjclfeo with mean 1 and finite 
variance. 

We will generally assume that the branching random walk is initiated by 
a single particle located at the origin £ Z. For a formal construction of 
a branching random walk following rules (A)-(B), see, for instance, [14]. 
The locations of the particles in generation n will be denoted by X„„ where 
i < Nn and N„ is the total number of particles in the nth generation. Our 
interest is in the distribution of the maximal displacement 

(5) M = maxM„ where Mn — max X„i. 

n>0 i=l,2,-,N„ 

Remark 1. Observe that only the particle locations at the end of each repro- 
duction step are taken into account in the definition of M: thus, for instance, 
if the initial particle at X04 = were to jump to site x = 1 and then produce 
no offspring, the maximal displacement would be M = 0, not M = 1. 

Many authors (e.g., | [T4| . IH) consider branching random walks in which 
the order of the reproduction and jump steps is opposite to that specified 
above. For our purposes it is more convenient to have the jump step 
precede the reproduction step, as this leads to more compact formulations 
of the main result and the nonlinear convolution equation ([TT]) below. It 
should be obvious that analogous results can be deduced for branching 
random walks of the type considered in |l4l, [6|, by conditioning on the 
first jump step. See Remark|2]in sec. l2.1l below. 

Theorem 1. Assume that the step distribution Frw - {fljtlfcez has mean 0, positive 
variance r^, and finite r-th moment for some r > 4. Assume also that the offspring 
distribution Few - {pk}k>o has mean 1, positive variance a^, and finite third 
moment. Then 

(6) P{M > x} as X ^ 00. 

This will be proved in section |2l The result can also be reformulated 
as a statement about the distribution of the maximal displacement of a 
branching random walk initiated by a large number n of particles at the 
origin. Such a branching random walk is just the superposition (sum) of 
n independent copies of the branching random walk in Theorem [H so the 
event that its maximal displacement is < y/nx is the intersection of the 
events that each of the n constituent branching random walks has maximal 
displacement < ^/nx. For fixed x > the target points yfnx ^ 00 as n — > 00, 
so the asymptotic formula ^ yields the following corollary. 
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Corollary 2. Let M" be the maximal displacement of a branching random walk 
initiated by n particles at the origin at time 0. Then for any x> Q, 



Theorem [T] is closely related to the main results of Kesten fT4], who 
considers critical branching random walk conditioned to survive for a large 
number of generations. (X. Zheng |19| has also considered the case where 
the random walk is assumed to have a "small" drift, but this leads to entirely 
different asymptotics.) Kesten shows that under the same hypotheses as in 
Theorem [ij for any fixed |S > 0, given that the branching process survives 
for j6n generations, the conditional distribution of maxj;<„ M/c/ yfn converges 
as n — > oo. He remarks on the result of Sawyer and Fleischman: 

We have not proved [equation in our setting. It is not clear 
at the moment how the methods of Sawyer and Fleischman, 
which rely on differential equations, can be carried over to 
the discrete setting; differential equations will have to be 
replaced by recurrence relations. 

Our main technical innovation will be to show how to relate these "recur- 
rence relations" to the differential equation (|3|. This will be accomplished 
by exploiting Feynman-Kac formulas. Our approach applies also to time- 
dependent Feynman-Kac problems, and leads in particular to information 
about the distribution of the random variable M„. As an illustration, we 
will prove in section |3] the following conditional limit theorem. 

Theorem 3. Under the hypotheses of Theorem [IJ the conditional distribution of 
Mn I ^fn, given that the branching process survives for n generations, converges 
weakly as n ^ oo to a nontrivial limit distribution G that depends only on the 
variances and of the offspring and step distributions. 

The scaling in this theorem is the same as that in the Dawson-Watanabe 
theorem (see, for instance, [10], ch. 1), which can be stated as follows. Sup- 
pose that n independent copies of the branching random walk are initiated 
at the origin £ Z. If particles are assigned mass 1 /n, and if time and space 
are re-scaled by factors 1/n and 1/ ^fn, respectively, then the corresponding 
measure- valued processes BRV\l„{t) converge weakly to super-Brownian mo- 
tion Xf. A similar theorem holds for the measure-valued processes BRW* (f) 
attached to branching random walk initiated by a single particle at the ori- 
gin, but conditioned to survive for n generations: under the same re-scaling 
of mass, time, and space as in the Dawson-Watanabe theorem, the measure- 
valued processes BRW*(f) converge weakly as n — > oo to a measure-valued 
process Yt- The law of this process Yt is related to that of the super-Brownian 



(7) 



n- 



lim P{M" > ^fnx} = 1 - exp{-C/x^} where C = -4-. 
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motion by the Poisson clustering representation: if N is a Poisson random vari- 
able with mean 1 and , Y^, . . . are independent copies of Yt then Y!i=\ 
is a super-Brownian motion. 

The weak convergence of the measure- valued processes BRW*(f) by it- 
self does not imply Theorem |3l because the location M„/ ^Jn is not a con- 
tinuous function (relative to the weak topology on measures) of BRW*(1). 
(Lalley l|T6l shows that in dimension 1 rescaled branching random walks 
converge to super-Brownian motion in a stronger topology than the weak 
topology implicit in the Dawson-Watanabe theorem. However, even this 
topology is too weak to make the normalized rightmost particle location a 
continuous functional.) Nevertheless, it is natural to wonder whether how 
the limit distribution G of Theorem |3] is related to the limiting measure- 
valued process Yt- The proof of Theorem |3] will establish that G is the 
distribution of the rightmost support point of the random measure Yi . 

Corollary 4. Under the hypotheses ofTheorem^(in particular, under the assump- 
tion that the step distribution Frw has finite rth moment for some r > A), 

(8) G(x) = P{Yi[x,c«) = 0}. 

Are 4-1- £• moments on the step distribution really necessary for the validity 
of our theorems (and Kesten's)? The following simple heuristic argument 
(which with a bit of work can be made rigorous) shows that 4 - e moments 
are not enough. Consider, for instance, the case where Frvv is the symmetric 
distribution on the nonzero integers with discrete density 

^^^^^^ = 2C(5-l)|xP-' 

which has infinite 4th moment. Conditional on the event that the branching 
random walk survives for at least n generations it will produce on the 
order of n^ particles. Each of these has conditional probability ~ Cln^~" of 
placing an offspring at distance n(i+'^)/2 to the right, where a - ell + ed - 46. 
Consequently, if a > then for large n the probability that all n^ particles 
are located in an interval [—A ^/n,A^/n] is vanishingly small. 



2. Maximal Displacement: Proof of Theorem [T] 

2.1. A Nonlinear Convolution Equation. For the remainder of the paper 
we shall assume that the offspring distribution few = {Pk}k>o arid the jump 
distribution Frw = {flxlxez satisfy the hypotheses of Theorem[ll in particular, 
f GW has mean 1, positive variance a^, and finite third moment, and Frw 
has mean 0, positive variance rf-, and finite 4 -I- £ moment. The maximal 
displacement M of the branching random walk is defined by equation 
and its (tail) distribution function will be denoted by 

(9) u{x) = P{M > x}. 
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Clearly, u{x) = 1 for all x < 0, and lirrij^oo u{x) = 0. 

We begin by showing that u satisfies a nonlinear convolution equation 
analogous to the Fleischman-Sawyer equation (|3]). This is obtained in the 
conventional manner, by conditioning on the first generation of the branch- 
ing process. Each particle of the first generation will give rise to its own 
descendant branching random walk, independent of its siblings (condi- 
tional on their locations), and in order that M < x it must be the case that 
the maximal displacements of all of the descendant BRWs, adjusted by their 
starting points, must be < x. This implies that for all x>l, 

oo 

(10) 1 - u{x) ^Yj^vYj p^^'^ ~ "(^ ~ y^^'- 

yeZ A:=l 

(Note that this is consistent with our convention regarding the sequencing 
of the dispersal and reproduction steps - see Remark [l] above.) Rewriting 
this equation in terms of u leads immediately to the following proposition. 

Proposition 5. u{x) satisfies the nonlinear convolution equation 

(11) u{x) = Y,'^kQ{u{x-k)l 

keZ. 

where 1 - Q(l -s) is the probability generating function of the offspring distribution 
Few, ihat is, 

CO 

(12) Q(s) Vii^ - for 0<s<l. 

i=0 

Remark 2. If the branching random walk used the alternative rule discussed 
in Remark[l](that is, particles first reproduce and then disperse) to construct 
the branching random walk, then equation l|TT]l would change as follows. 
Writing M for the maximal displacement of this branching random walk, 
and u{x) = P{M > x}, we would have 

Hx) = Q(Yj a^uix - k)). 

/ceZ 

Comparing this with equation ((TT]) . we see that 

u{x) = Q(u{x)). 

Since the Taylor expansion of Q is Q(s) = s - o^s^/2 + 0{s^), it follows that 
u(x) and u(x) go to as x — > oo at the same rate. 

2.2. A Discrete Feynman-Kac Formula. Our goal now is to analyze the 
asymptotic behavior of solutions to the nonlinear convolution equation l|TT]) 
as X — > oo. To accomplish this, we will show that solutions of (|TT]) can 
be represented by formulas of "Feynman-Kac" type. Henceforth, we shall 
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denote by W„ the random walk on TL whose step distribution is the reflection 
of the step distribution f rw in the underlying BRW , that is, 

(13) P(W„+i - W„ = y I W„_i, ...)=a.y. 

We shall use superscripts and £^ to denote the initial point Wq = x of the 
random walk W„. 

Define 

(14) }i{s) = s - Q(s) = o^s^jl + 0(s3) and 

H(s) = /j(s)/s = 0^5/2 + 0(s2). 

It is easily checked that H(s) is increasing for s € [0, 1], and satisfies H(0) = 
and H(l) = po- 

Proposition 6. Under P^, the process 

n 

(15) y„ = ([](l - H(M(Wy)))) • M(W„) 

7=1 

/s a hounded martingale with respect to the natural filtration generated by the 
random walk {W„}. Here we use the convention that the empty product n°=i 
equal to 1. 

Proof. The random variable {Y,,} is clearly bounded: < < 1. To prove 
that the sequence is a martingale we appeal to the fundamental convolution 
equation (|TT]l. which can be rewritten in terms of the function h as 

/ ^ 

(16) ^ fl)tu(x -k) - u{x) - ^ aj^h{u{x - k)). 
^fceZ ] keZ 

Using this, we compute 





1 mu) 




n 
7=1 


HiuiWj)))) 


E'((l - H(m(W„+i))) • M(W„+i) 1 {Wj]1^^ 


n 

7=1 


HiuiWj)))) 


£-^((w(w„+i) - h{u{w„^,))) 1 { w^};-^,) 


n 

7=1 


H{u{Wj)))) 


fl,w(W„ - fc) - ^ «)c/i("(W„ - k))) 
fceZ teZ 


n 
7=1 


H{uiWj)))) 


U(W„) = Yn, 



where the second to last equality uses the fundamental equation ((T6|) . □ 
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Corollary 7. For each y eZ, define 

(17) Ty = mm{n > | W„ < y}. 

Then for all x, y £ Z, 



(18) u{x) = E 



y 

■X 



]\[l-H{u{]Nj))) 



"(W.„). 



Proof. Since the random walk W„ is driftless it must be recurrent, and hence 
Ty is finite. Since the martingale Y„ of Proposition [6] is bounded, Doob's 
optional sampling identity applies, yielding ((Ts)) . Note that if y < then 
M(W^y) = 1, since W^^ < 0. □ 



2.3. Scaling Limits. Using the Feynman-Kac representation ([TSt we will 
show that the function u, properly re-normalized, converges to a function 
that satisfies the Fleischman-Sawyer equation (|3]). Because we do not know 
a priori that the function u has a proper scaling limit we must work with 
subsequential limits. However, since u is monotone and satisfies < u < 1 
it follows that for any y > there exist sequences X)t — > oo such that 

no\ A.( \ V "fa + y/ V»fa)) 

^(^)^=,Sl^ ^fe) 

exists. Clearly, the limit must satisfy < (/)(y) < 1, and if y = the limit is 
(f)(fS) — 1. By Cantor 'sdiagonalization argument, any such sequence X)t must 
have a subsequence, which we also denote by X)t, such that the convergence 
(fT9l) holds for all rational y > 0. 

Proposition 8. For any sequence — > oo suc?z f/iflt ([T9l > /lo/ds /or aZZ rational 
y > 0, t?je limit function (p{y) extends to a continuous, non-increasing, positive 
function of y e [0, oo). Hence, the convergence ([T9l > fzoWs uniformly for y in any 
compact interval {Q,A\. 

Proof. Fix < yi < y2, both rational, and for ease of notation write z, = 
y,7 ^u{x) for / = 1, 2. Fix a sequence X/^ — > oo along which ([191) holds for all 
rational y. To avoid a proliferation of subscripts, we shall write lim^^oo to 
mean convergence along the subsequence x^. By Proposition[6]and Doob's 
optional sampling theorem, 

U{X + Z2) 

(piyi) = lim 



u{x) 



= lim E^+"2 

1'— >oo 
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where Zz - t(z) = mini/ > | Wy < z}. Using the expansion H{u) ~ as 
M — > we obtain 



(p{y2) = lim 



= limE^+^2-"i 

X—foo 



u{x) 



[](l-H(u(W,+zi))) 



> limE^+^2-"i |exp|-yT(x)u(x)| 



u{x + Zi) 



m(x) 

The last inequality holds because Wj > x for all ; < and < x. 

By the invariance principle (here we use the assumption that the step 
distribution of the random walk W„ has mean and finite variance), as 

X — > CO, 

D{TrU{x) I Wo = X + Z2 - zi) = D{tou{x) | Wq = Z2 - Zi) 

^2)(T««|Bo = y2-yi) 

where Bf is a standard Brownian motion started at 1/2 ~ J/i and is the 
first hitting time of by Bf . Hence, for any e > 0, when 1/2 - yi is sufficiently 
small and x is sufficiently large, exp{-a^TxM(x)/2} > 1 - e with probability 
1 - e. Consequently, 



limE^+"2-^i (exp(-^T^M(^)) 

X^oo \ ^ 2 ' 

>(l-e)^lim"(" + ^i^ 



u{x + Zi) 
m(x) 



u{x) 



= (l-e)2(/,(yi). 



This proves that Ve > 0, if Z2 - Zi is sufficiently close to zero then (p{yi) > 
(piyi) ^ (1 - £)^(/^(yi)- Therefore, (p is continuous. □ 

Our aim now is to show that there is only one possible subsequential limit 
function (p, and that it satisfies the Fleischman-Sawyer differential equation. 
To accomplish this we will use the discrete Feynman-Kac formula (|l8]) and 
the invariance principle to show that any subsequential limit (p satisfies the 
following Feynman-Kac formula. 

Proposition 9. Assume that the step distribution {fljtltez of the random walk {W„} 
has finite r-th moment for some r > 4. Then any subsequential limit (p{y) specified 
by ((19)) satisfies 



(20) (Piy) 



2 />T™ 

= EW')exp(-y (P{i]Bt)dt), 
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where under P^'l the process Bt is a standard Brownian motion started at Bq - y/rj 
and is the first hitting time ofO by Bf. 

The need for finite A + e moment stems from the fact that in general 
the random walk W„ will overshoot at the first passage time t{0). The 
renewal theorem (cf. [llj or L18J) implies that as the initial point Wq = x — > 
oo, the distribution of the overshoot W^(o) converges weakly provided the 
step distribution of the random walk has mean zero and finite variance. 
The number of finite moments of the limiting overshoot distribution is 
determined by the number of moments of the step distribution, as follows. 

Lemma 10. If the step distribution of {W„} has finite r-th moment, then the 
limiting overshoot distribution has finite (r - l)-th moment. 

Proof. Consider the ladder variables 

Ti = min{n > | W„ < Wq}, Zi = Wj, - Wq, 

T2 = min{n > Ti | W„ < Wj^], Zi = Wt^ - Wt,, 

and so on. Clearly, the first passage time t(0) must be one of the ladder 
times Tj. Moreover, the ladder steps Z,„+i - Zm are i.i.d., and by exercise 
6, p. 232 of IflSl , the random variable Zi has finite absolute rth moment 
if the step distribution Frvv has finite absolute (r + l)th moment. The key 
observation is that for any fl < and any x > 1, 

X 

P^(W,(o) <a)^Yj G{x}k)PHZi <a-k} 
k=l 

X 

<Y,P'iZi<a-k] 
k=l 

X 

=2]^p°{Zi<fl-fc} 

k=l 

where G(x; k) is the probability under P^ that the random walk W„ will visit 
the site k at one of the ladder times T,. Thus, if S has the limiting overshoot 
distribution, then for any fl < 0, 

oo 

P{S <a)<J^ P{Zi <a-k). 

k=l 

This inequality together with the earlier observation about moments of the 
ladder variable Zi implies that if E'^|Wi|'' < oo then S has at least r - 2 
moments. □ 

Lemma 11. If the step distribution f rw satisfies the hypotheses of Theorem^then 
along any sequence x = Xjt — > oo such that ((191) holds uniformly on compact sets, 

uiW^^ + Xk) \ ^ ^ 

u{Xk) I 



(21) 



limEf/' 
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Proof. As in the proof of Proposition [H] we will omit the subscript k on x^^ 
and write limj^oo to mean convergence along the subsequence xi^. We also 
write z-yl ^u{x). By Proposition|8l 

u{x + z) 
lim iim — ; — = 1, 

which implies that for any a £ (0, 1 /2), 

u(x + u(x)~^'^'^) 
km — = 1. 

By the monotonicity of u, 

\ u{x) I 

where 

A = A{x) = {W,(o) > -w(x)-i/2+«}. 
By Chebyshev's inequality, for any r > 2, 

P'{A') < E|W,(o)|'-2m(x)(''-2)(5-«), 

and by Lemma [lOl if the step distribution Frw has finite rth moment then 

E|^t(0)I'~^ < Consequently, 

11 < C{uix)f2-")('-2)-\ 

Since by hypothesis the step distribution Frw has finite rth moment for 
some r > 4, the constant a > can be chosen so that the exponent in the last 
displayed inequality is positive, and so it follows that quantity 11 converges 
to as X — > oo. On the other hand, lim^^oo P^{A) - 1, and on the event A 
the integrand in quantity I is bounded by 

m(x-m(x)-i/^+«) _^ ^ 
u{x) 



and so quantity I converges to 1 as x — > oo. 



□ 
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Proof of Proposition\9\ Once again write z - y/ ^u{x). According to Proposi- 
tion[6]and the Optional Stopping Theorem, 

u(x + yl y/u(x)) 
u{x) 

= E«n(l-H(.,(W,)))^ 



To 



w(W^o + x) 



(22) 



|"(Wxo+X) 



= exp(2^ log(l - H{u{Wj + x)))-^ 

;=1 ^ ^ 

= E^exp(2^(--«(W, + X) + 0{u{xf))] 

= exp|-- 2^ u(W;- + X) + ToO(u(x)^)) J° 

;=1 ^ ^ 

The error term 0(w(x)^) is bounded in magnitude by Cu{xf- for some finite 
constant C not depending on x, by virtue of our standing hypothesis that 
the offspring distribution few has finite third moment, which ensures that 
H{u) has finite third derivative at u = 0. 

The invariance principle implies that as x — > oo the distribution of the 
process V"W^t/M(.T)/'] under converges weakly to that of a standard 
Brownian motion Bf started at Bq = y/U- Consequently, the distribution 
of the renormalized first passage time m(x)to converges weakly to that of 
Tq^, and hence the error term toO(u(x)^) converges in distribution to 0. 
Moreover, along any sequence x = x/^ ^ oo such that the convergence ([l9l ) 
holds uniformly for y in compact intervals. 

To To ^t:BM 

y w(Wy + x) = m(x) y M(Wy + x)/u(x) ^ I (p{riBt)dt. 

^0 

Therefore, by Lemma [TTl 

lim expi-- y u{Wj + x,) + ToO(u(Xfc)^) ;° 
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This together with the convergence ((19)) and the chain of equalities ((22]) 
proves that (p must satisfy the Feynman-Kac formula ((20|) . □ 

Corollary 12. Under the hypotheses of Proposition^ 

u(x + yl Ju{x)) (ay \^ 
(23) lim " + ~^{y). 



\ V67] 



Proof. It suffices to prove that there is only one possible subsequential limit 
function ( (191) , and that this limit is the solution of the differential equation 

(24) (\>"{y) = a^iyfln^ 

that satisfies (p{Q) - 1 and limy^oo (p{y) - 0. This follows from the Feynman- 
Kac representation ((20l ) of subsequential limits and Kac's theorem (cf., for 
instance. Theorem 3.7, ch. II of [5]), which implies that for any positive, 
bounded, continuous function 1/ : [0, oo) — > R the function 



V.(y) = £J'/''exp V{^Bt)dt 




is the unique bounded solution of the differential equation ip" - o^Vip/r]^ 
satisfying ip(0) = 1. □ 

2.4. Proof of Theorem (H To complete the proof of Theorem [T] we must 
show that 

(25) lim w{x) = := Sif/o^ where iv{x) := x^u{x). 



This we will deduce from the asymptotic scaling law (|23l >, which in terms 
of the fimction w may be rewritten as 

(1+iSy)' io{x{l + y/^/Mx))) ^ 
(26) lim —= — = 1. 



Lemma 13. limsupK;(x) < oo. 

Proof. First we show that lim inf^^oo w{x) < oo . If not, then for all sufficiently 
large k the function w would assume the value k for a last time, atZk :- sup{x | 
io{x) < k} < oo, and zv{Zk) - khy the continuity of w. Since w{z) > k for all 
z > Z|f, it would then follow that for any y > 0, 

But 



> 1 



(1 + Byf , 
lim = (1 + |3y)2 > 1, 



so this would contradict the relation ((261). 
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Next, suppose that limsup^^^^ zv{x) - oo. Since liminfj^oo ^v{x) < oo, by 
the preceding argument, for any sufficiently large A there exist sequences 
Xk> Xk ^ CO such that w{Xk) - A^/fi'^ for all k but lim^t^co u){xk) - oo. Define 
ykhy 

Xk = Xkil + i3i/)tM); 

then 

But by construction, 

wixkil + yk/ V^fe))) 

= > oo, 

so again we have a contradiction to (|26l l. □ 

Proof of Theorem^ We must establish that lim^^oo W = l/jS"^- Let x;^ ^ co 
be any sequence along which iv{xk) converges, and denote L = limj-^oo w{xk). 
By Lemma [l3l the limit point L must be finite. Fix y > 0, and choose 
- < such that 

Xk = Xkil + y/ ^|w{xk)) = Xk + y/ ^|u{xk). 

Such points Xk > must exist, provided Xk > y, because the right hand 
side increases continuously with Xk on the interval [0, Moreover, for 
any fixed y > the sequence Xk{y) — > oo, since X)t — > oo. By going to 
a subsequence of Xk if necessary we may arrange that the sequence w{xk) 
converges, with limit K - K{y) < oo. The fundamental relation ((26l) (using 
X - Xk) now implies that 

(l + |gy)^ L _ {^^L + p^^Lyf _^ 
(l + y/VK)2 (VK + y)2 

and hence 

(27) VL + (|3VL-l)y. 

By Lemma[l3l the set of limit points K(y) is bounded above. Since this must 
hold for every y > 0, it follows from ((27|) that jS^L < 1. On the other hand, 
for every y £ K every limit point K{y) must be nonnegative, so (|27)l implies 
that jS^L > 1. 

□ 

3. Conditional Limit Theorem 

3.1. Space-time Feynman-Kac formula. Assume throughout this section 
that Mn is the rightmost particle location in the nth generation of a branch- 
ing random walk satisfying the hypotheses of Theorem[Tl The unconditional 
distribution function of the random variable M„ satisfies a time-dependent 
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nonlinear convolution equation similar to the time-independent equation 
([Tl]) satisfied by the distribution of the maximal displacement random vari- 
able M. Specifically, if 

(28) Vn{x) = v{n,x):=P{Mn>x}, 
then for every n >1, 

(29) Vn{x) = ^ akQ{Vn-iix - y)), 

where 1 - Q(l - s) is the probability generating function of the offspring 
distribution (see equation ((12)) ). The objective of this section is to analyze 
the asymptotic behavior of v, and in particular to show that nv{n, [x ^Jn\) 
converges as n — > oo, for any :t £ R, to a distribution function that de- 
pends only on the variances of the offspring and step distributions of the 
branching random walk. The strategy will once again be to represent the 
solution of (|29l) by a discrete Feynman-Kac formula, and then to show that 
after an appropriate rescaling the Feynman-Kac expectations converge to 
a Feynman-Kac expectations for Brownian motion. As in sec. |2l denote by 
W„ a random walk with step distribution dlSl) , and indicate by P'^ the law 
of the random walk with initial point Wq = x. Then essentially the same 
arguments as in the time-independent case prove the following assertion. 

Proposition 14. For each n the process 

k-l 

(30) Zf ) = v„.k{Wk) 11(1 - H{Vn.j{Wj))) 
is a martingale under P^. Consequently, for any n <m, 

m-j 

(31) V„r{x) = E^Vn{Wm-n) - H{v„,-j{Wj))). 

3.2. Monotonicity, tightness, and scaling limits. For each n > the func- 
tion Vn{x) is non-increasing in x, with limit as a: — > oo and limit z;„(-oo) = 
P{C > n} as X — > -oo, where C is the extinction time of the branching random 
walk. By Kolmogorov's theorem on the lifetime of a critical Galton- Watson 
process (see IH, ch. 1), as n ^ oo, 

(32) Vn{-^) = P{C > n} ~ 
Hence, the functions nv„{x) are uniformly bounded. 

Lemma 15. Under the hypotheses of Theorem^ the family ofrescaled distribution 
functions nv„{x yfn) is tight, that is, 

(33) lim sup nP{Mn > x ^fn} - and lim sup nP{M„ < x ^fn} - 0. 

X^oo ^^^^ x^-oo 
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Proof. The first of these follows directly from Theorem [H because M„ < M. 
The second follows from Theorem[l]by reflection of the branching random 
walk in the origin. □ 

Remark 3. The hypothesis that the step distribution Frw of the branching 
random walk has finite 4 + £ moment is used here in an essential way. If 
Frw has infinite 4 - f moment for some e > then Lemma Hsl need not be 
true: see the discussion at the end of sec.|2l 

Lemma 16. For any £ > there exists 6 > such that if\x-y\<6 yfn and 
n <m < n{\ + 6) then 

(34) ^ - 1 < 



v„{x) 

Proof. This follows routinely from the Feynman-Kac formula dSTl) and Donsker 's 
invariance principle, since sup^y„(x) — > as n — > oo and H{u) ~ a^ujl as 
u ^ 0. □ 

Corollary 17. Any sequence of positive integers has a subsequence nj^ ^ oo along 
which the functions {t,x) i-> nv[nt]{x^fn) converge uniformly for 1 < t < A and 
x € [-00, (xi],for any A < oo. The set of possible limit functions (p{t, x) is compact 
in C{[1,A] X R, and for each t the function (p{t,x) is non-increasing in x, with 

(35) lim(p(f,x) = and lim t<p{t,x) = 2/o'^. 

X — >oo X — >— oo 

Proof. All of the assertions except the limits ((35)) follow from Lemma[T6l The 
limit relations i35\) follow from the tightness of the family nvn{x ^/n). □ 

Corollary 18. Let (p{t, x) be any subsequential limit of the functions {t, x) i-> 
nv[nt]{x yln),for 1 <t and x e [-oo, oo]. Then cp satisfies the identity 

(36) (p{t + l,x)^ E^I\{1, r]Bt) exp !^-^ (p{t + 1 - s, 7]B,) rfsj , 

where under py the process Bt is a standard Brownian motion started at y. Conse- 
quently, (p satisfies the partial differential equation 

(37) ^ = — ^-aV for t>l andxeR. 
ot 2 dx^ 

Proof. The integral representation (|36l ) follows from the discrete Feynman- 
Kac formula jST] ) and the invariance principle by virtually the same argu- 
ment as in the proof of Proposition |9l The differential equation dSTl) follows 
from the integral representation by Kac's theorem. □ 

It should be noted that the only a priori bounds on the fimctions Vn 
needed to deduce the existence of subsequential limits (pit,x) are those 
in Lemmas [T5l - [T6l and Corollary [iTl These use only the crude estimate 
M„ < M :- max„>iM„ and the results concerning the tail behavior of the 
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distribution of M proved in section 12.21 To prove that there is only one 
possible subsequential limit function (p{t,x) we will need the following 
stronger a priori bounds on the functions y„. 

Lemma 19. 

(38) lim supnx^y„(x Vm) = and 

n>l 

(39) lim inf nx^{l - v„{x ^fn)) = 0. 

x->-oo n>l 

The proof is deferred to section |331 below. 

3.3. Super-Brownian motion and solutions of (|37)) . At first sight it might 
appear that Corollary [18] sheds no light at aU on the question of uniqueness 
of scaling limits, because in general the solution to the partial differential 
equation (|37)) will depend on the initial condition (p{l, x). However, we will 
show, using the a priori bounds in Lemma[T9l that in fact CoroUaryllSlimplies 
that there can be only one scaling limit, and thereby complete the proof of 
Theorem m The key is the fact that solutions of (|37)) determine - and are 
determined by - the law of super-Brownian motion, by the following duality 
formula. For ease of exposition, assume henceforth that rf- = 1. (There is no 
loss of generality in this, because solutions of (|37)) can be rescaled.) 



Proposition 20. Let X^" be a super-Brownian motion with branching parameter 

and initial mass distribution X^'' = 6x, a unit point mass at location x £ IR. 
Let cp{t + \,x) be the solution of the evolution equation (|37l ) with initial condition 
(p{l, x) - q){x). Then 

(40) (pit + 1, X) = - log E exp{-(Xf% (p)}. 

Here (p, (p) denotes the integral ofcp against the measure p. 

Proof See, for instance, [lOl, Cor. 1.25. □ 



Exploitation of the duality formula (|40|) will require some elementary 
properties of super-Brownian motion. First, super-Brownian motion is 
equivariant under spatial translation: xj** can be obtained from X^° by 
translating each of the random measures " by x to the right. Second, super- 
Brownian motion satisfies a scaling law: if Xf is a super-Brownian motion 
with initial mass distribution Xq = ^^-^ ^^^^^ the rescaled measure- valued 
process Xf defined by 

<Xf,/) = <A-iXAf,/(VA-)) 

is a super-Brownian motion with initial mass distribution Xq = bx- Third, 
super-Brownian motion is infinitely divisible, in the following sense: the 
superposition of m independent super-Brownian motions with initial mass 
distributions v, is a super-Brownian motion with initial mass distribution 
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TjILi '^i- This implies that the total mass |Xf | at time t evolves as Feller diffusion, 
and so P{Xt 0} ~ C/t as t — > oo with C = where is the branching 
parameter (cf. [?], sec. 2.1, Theorem 1 and Example (iii)). It follows, by 
(|40)) . that any solution (p{t + 1, x) of ((37|) whose initial condition (p(x) - (p{l, x) 
is non-increasing in x and satisfies the boundary conditions (|35]l . Finally, 
the scaling and infinite divisibility properties imply that super-Brownian 
motion conditioned to live for time at least t has a limit law: if Xf is a 
super-Brownian motion with initial state Xq = 6o then for any bounded, 
continuous test function / : R — > IR, 

(41) lim D{t-\Xt, /( VI •)) I ^ 0) = D«y^, /)) 

where Yt — is the super-process specified in Corollary HI 

Construction of the limit process. The existence of the weak limit (|4T]l is 
implicit in the "Poisson cluster" representation of super-Brownian motion 
(cf. [10], ch. 1) but in fact weak convergence will not by itself suffice for the 
arguments to follow, so we now sketch a construction of super-Brownian 
motion in which the random measure arises naturally as an almost sure 
limit. First, observe that infinite divisibility implies that a super-Brownian 
motion Xf with initial state Xq = 6x can be decomposed as 

Xf = X[ -h X", 

where X^ and X" are independent super-Brownian motions started from ini- 
tial measures X^ = Xg = 6x/2. This decomposition process can be iterated, 
so by standard arguments there exist, on some probability space, countably 
many independent super-Brownian motions X"'"', where m = 0, 1, 2, . . . and 
n = 1,2, ■■ ■ , 2'", such that for every pair (n, m) the process X"''" has initial 
state X^''" = 6J2'" and 

By the scaling property, the probability that any one of the mth generation 
processes {X"''"}f>o survives to time 1 is the same as the probability that the 

Oth-generation super-Brownian motion {X^f'^]t>Q survives to time 2"', which 
is ~ C/2'", where C = 21 o'^. Hence, if F„, = {n \ X"'"" + 0} then |F,„| converges 
in law to the Poisson distribution with mean C. By construction, the random 
variables |F,„| are nondecreasing in m, so it follows that in fact the sequence 
{|Fm|}m>o is eventually constant, with Poisson limit N. By re-indexing the 
super-Brownian motions in each generation m, we can arrange that for all 
sufficiently large m the processes {X"'"'}t>o satisfy 

X^' = XY""-^ for all 1 < n < N, 
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and only those processes {X^"''"^}t>o for which n < N survive to time t - 1. 
Conditional on the value of N, each of the random measures 

X"'°° := lim X"''" 
is an independent version of Yy 

□ 



3.4. Asymptotic behavior of scaling limits. 

Corollary 21. Let (p(t, x) be any subsequential limit of the functions {t, x) i-> 
nvi„tj{x yln),for 1 <t and x e [-oo, oo]. Then for any xeW., 

(42) lim 2"'(p(2'" + 1, xT'l^) = ^P{Y^J-oo, 0] ^ 0} := 2G{x)/o^. 

furthermore, for each x £ R f/ze convergence ((42)) holds uniformly over the set of 
all possible subsequential limits (p{t, x). 

Proof. Write y = ^/ix, and abbreviate (p{l,x) = (p{x). Since P{X^'' ^ 0} ~ 
2/(a^0 as f — > oo, the duality formula ( |40l) implies that 

t(p{t + 1, y) = -t log E exp{-<xf^ <p)} 
~fE(l-exp{-<xf",<p)}) 
= fE(l-exp{-<xf",<p)})l{xf"^0} 
- (2/a2)E((l - exp{-<xf^ <p)}) | x'" t 0). 



The last expectation can be rewritten using the scaling property of super- 
Brownian motion: 

E{{l-exp{-<x'",<p)})|x'"^0) 

(p(z)d(x'^(z))}|xf" ^0) 

OO 

(pi Viz) d(rixf"( Vfz))} I x^^ ^ 0) 

oo 

Xoo 
<p(Viz)d(x;'^^*(z))}ix;'^/'^o). 

By the construction sketched in the preceding subsection, versions of the 
super-Brownian motions X^'^*, for t - 2'", conditioned to survive to time 
t - 2"' can be constructed on a common probability space along with a 
version of the random measure in such a way that X^""^^ = Y^ for all 
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large m. Hence, as ^ — > oo through powers of 2, 

oo 

(p(Vfz)d(y^(z))}. 

oo 

The result now follows from Lemma [191 and the dominated convergence 
theorem. To see this, observe that the exponential in the last expectation is 
bounded above by 1, because the function cp is nonnegative. By Lemma [191 
for each z > 



lim t(p{ yftz) - and 

t — >oo 

lim t<p{- yftz) - -oo 



f — >OQ 



uniformly over the set of all possible subsequential limits (p{y) = lim nvn{ "0^3/) • 
Therefore, by dominated convergence, as i — > oo through powers of 2, 

Xoo 
(p{ Vtz)d(y?(z))} = P{supp(Y^) c (0,oo)}. 
oo 



□ 



3.5. Proofs of Theorem[3land CorollaryHl It suffices to show that 

2 

lim nv„{xyfn) - V{x) :- -^G{x), 

n— >O0 

where G(x) is defined by (|8ll. By Corollary [I7l subsequential limits exist, 
and by Corollary[l8lsubsequential limits must satisfy the partial differential 
equation ([37l) . What must be shown is that the only possible limit is the 
function V. 

Suppose then that there is a subsequence rij^ ^ oo along which nv„{x) — > 
U{x) for some function U. By CorollaryllTlthe sequence n/^ has a subsequence 
Uj — > oo along which the functions {t, x) i-> nv[„t] {x ^fn) converge uniformly 
for \ < t < A and x e [-oo, oo], for any A < oo. By rescaling time, we can 
extract yet another subsequence n, along which the the functions (f, x) i-> 
V^) converge uniformly for t € [2~'",2'"] and x £ [-00,00], for any 
m >1. Denote the limit function by <p{t, x). 

By construction, <p{l,x) = U(x)} moreover, each section cp(2~'",x) is a sub- 
sequential limit of the functions x i-> 2~"'ny2-"')i(^ V^2"'/'^). Consequently, 
Corollary [211 implies that U - V. 



□ 
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3.6. Proof of Lemma [191 The proof will use Theorem [T] and the Dawson- 
Watanabe theorem. Denote by £,],£,^,... the (counting) measure-valued pro- 
cesses associated with independent branching random walks each started 
by a single particle at the origin at time 0, and set 

n 

(43) = 2^ 

i=l 

Thus, the process {S"}t>o is a branching random walk initiated by n particles 
all located at the origin. (Here we view the branching random walks as 
continuous-time processes that are constant on time intervals [m,m + 1), 
with jumps at integer times m.) The Dawson-Watanabe theorem asserts 
that the process {S"}t>o, after rescaling, converges in law as n — > oo to super- 
Brownian motion Xf with initial mass distribution Xq = 6o- In particular, 

(44) h"J^/^■)^Xt 

where the weak convergence is in the Skorohod topology on the space of 
cadlag measure- valued processes. The limiting super-Brownian motion has 
local branching rate t]^ and diffusion coefficient a^. 

Super-Brownian motion Xt in one dimension has the property that with 
probability one, for each f > the random measure Xf is absolutely contin- 
uous relative to Lebesgue measure, with jointly continuous density X{t,x). 
The density X(f, x) is jointly continuous except at f = 0, and for each t > 
has compact support in a:; as f ^ the support contracts to the point 0. 
Super-Brownian motion dies out in finite time, so 

:= sup{x £ R : I Xt[x, oo) dt > 0} 
Jo 

is well-defined, measurable, and finite. Denote by Tx the infimal time that 
Xf[x, oo) > 0; then the event > x coincides a.s. with < oo. The path- 
continuity properties of the density X{t, x) imply that for any e > and any 
compact interval [xi,X2] not containing there exist < 6 < A < co such 
that for all x e [xi, X2], 

(45) P(6 < < A I T.^ < 00) > 1 - £. 
Proposition 22. For any x > 0, 

(46) P{M^ > x} = P{tx < 00} = 1 - expl-C/x^} where C = 

Proof. By a theorem of Dynkin (see, e.g., ITOl , Ch. 8) the function 

u{x) = - log P{M^ < x} 

is the unique solution of the differential equation u" — ri^v?- ja^ with bound- 
ary conditions m(0) = co and 1/(00) = 0. □ 
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The crucial point of Proposition |22] is that the distribution (f46|) coincides 
with the limit distribution in Corollary |2l As noted earlier, the weak con- 
vergence (f44| does not by itself imply that the maximal displacement M" 
of the branching random walk S" converges weakly (after rescaling) to M^, 
because in the Dawson-Watanabe scaling (f44|) individual particles receive 
vanishingly small mass as n — > cx), but it does imply that 

liminf P{M" > y/nx} > P{M^ > x}. 

Thus, Proposition |22] and Corollary |2] together imply that the event {M" > 
'\Jnx} is almost entirely accounted for by sample evolutions in which a 
large number (order n) of particles reach yjnx. Furthermore, by (f45|) , they 
imply that for large n the event {M" > ^Jnx] is mostly composed of sample 
evolutions in which particles of the branching random walk reach [ ^fnx, oo) 
during the time interval [n6, nA], where < 6 < A < oo are as in (|45l >. 
Following is a formal statement of this observation. 

Corollary 23. Denote by M" the location of the rightmost particle of the branching 
random walk S" at time t. Then for any £ > and any compact K c (0, oo) there 
exist constants < 6 < A < oo such that for all x e K and all sufficiently large n, 

(47) nx^P{Mf„ > yfnx for some f £ [6, A]} > (1 - £)(1 - expl-C/x^}). 

Proof of Lemma^9\ Recall that v„{y) is the probability that the maximal dis- 
placement of the branching random walk E,^ exceeds y. Since the reflection 
of this branching random walk is again a driftless, critical branching ran- 
dom walk, the inequality dSS)) implies inequality (l39l ), so it suffices to prove 
([Ml/ that is, 

(48) lim sup nx^P{^l[ yfnx, oo) > 1} = 

^^'^ n>l 

We will accomplish this by contradiction. Suppose that there exist y > 
and sequences nj^, — > oo along which the left side of (|48l) remains bounded 
below by y. Let m = mj^ = [nj^x^] and 9 - 9k = 1/^^/ and consider the 
branching random walk S"' initiated by m particles at the origin; then our 
hypothesis would imply 

(49) P{MZe > V^} > y' = ^(1 - exp{-C/y}) > 

along the subsequence m = m^. We will show that ( [49l) is impossible. 

Denote by the event that M'^q > ^/m. We begin by observing that on 
this event the total number of particles in the interval [|S ^/m, co), for any fixed 
|3 > 0, must be small relative to m. This follows from the Dawson-Watanabe 
theorem: since 9 = 9^ ^ 0, for any a > the chance that the limiting 
super-Brownian motion puts any mass on the interval [jS, co) at some time 
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t < 6 is vanishingly small, and so for any a, e > 0, ii k is sufficiently large 
then 

(50) P{SZe[^^,^)>am} <e. 
By reflection, it foUows that 

(51) P{SZg{-^,-liyI^]>am}<e. 

Now fix £ > small, let K = [1/2,2], and let < 6 < A < oo be as in 
Corollary |23l Since 6k = x^^ ^ as k ^ oo, eventually Ok < 6/2. By 
((5T|) . with probability at least 1 - 2e there will be fewer than am particles 
outside the interval [- y/m/2, y/m/2]. Consider the component branching 
random walks initiated by these particles at time Om: if a <^ 6/2 then by 
Kolmogorov's theorem on the extinction time of a critical Galton- Watson 
process, the (conditional) probability (given the history of the branching 
random walk up to generation [6m]) that one of these branching random 
walks survives to time 6m is vanishingly small (as a — > 0); in particular, for 
suitable a this probability will be < £. Thus, except with probability not 
exceeding 3e, on the event A"^' the only particles that will survive to time 
6m are descendants of particles located in [- ^/m/2, ^/m/2] at time 6m. 

Next, consider the total number of particles in the branching random 
walk at time dm. By Feller's theorem, the processes N'//^^/m converge in law 
to a Feller diffusion Ft started at Fq = 1. Since diffusion processes have 
continuous paths, and since 6 - 6^ ^ 0, it follows that with probability 
approaching 1, 




as A: — > oo. In particular, the probability that NJJ^g > (1 + e)m will eventually 
be smaller than e. 

Finally, consider the event F^^ that > Vmfor sometime f € [6, A]. For 
F^^^ to occur, one of the particles alive at time md must have a descendant 
that reaches the interval [ ^/m, co) after time 6m. The probability that such a 
particle is not located in [- ^/m/2, ^/m/2] is less than 3e. Moreover, in order 
that a particle located in [- ^/m/2], ^/m/2 at time 6m has a descendant in 
[ ^/m, oo), the intervening trajectory must travel a distance at least ^/m/2, 
and by Theorem [T] the (asymptotic) probability of this is less than ^C/m, 
where C = 6rf- /o^. But the probability that the number of particles in 
[- yfm/2\ Vm/2 at time 6m exceeds (1 + e)m is less than e, and hence, for 
large k, 

P{Fl^\A'^) < 4£ + 1 - (1 - 4C/?«)<^+'>'" « 46- + 1 - exp{-4C(l + £)}• 
It follows that 

P{{FlJf^A^)>{l-e-^^)P{A^,). 



24 



STEVEN P. LALLEY AND YUAN SHAO 



Since the event IS/T" > s/m contains F^'^ U A'^, we now have, by Corollary |2l 

P{M'" > Vm} ~ 1 - e"^ 

> (1 - £)(1 - e-^) + (1 - e-^^)y' 

Since £ > can be chosen arbitrarily small, this is impossible, so we have 
arrived at a contradiction. 

□ 
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